Super-radiance, Berry phase, Photon phase diffusion and Number squeezed state in 

the U(l) Dicke ( Tavis-Cummings ) model 
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Recently, strong coupling regimes of superconducting qubits or quantum dots inside a micro- wave 
circuit cavity and BEC atoms inside an optical cavity were achieved experimentally. The strong 
coupling regimes in these systems were described by the Dicke model. Here, we solve the Dicke 
model by a l/N expansion. In the normal state, we find a yN behavior of the collective Rabi 
splitting. In the superradiant phase, we identify an important Berry phase term which has dramatic 
effects on both the ground state and the excitation spectra of the strongly interacting system. The 
single photon excitation spectrum has a low energy quantum phase diffusion mode in imaginary 
time with a large spectral weight and also a high energy optical mode with a low spectral weight. 
The photons are in a number squeezed state which may have wide applications in high sensitive 
measurements and quantum information processing. Comparisons with exact diagonization studies 
are made. Possible experimental schemes to realize the superradiant phase are briefly discussed. 



Recently, several experiments Qj successfully achieved 
the strong coupling of a BEC of N ~ 10 5 87 Kb atoms 
to the photons inside an ultrahigh-finesse optical cav- 
ity. In parallel, strong coupling regime was also achieved 
with artificial atoms such as superconducting qubits in- 
side micro-wave circuit cavity and quantum dots inside 
a semi-conductor micro-cavity system [3|]. In these ex- 
periments, the individual maximum coupling strength g 
between the ( artificial ) atoms and field is larger than 
the spontaneous decay rate of the upper state 7 and the 
intra-cavity field decay rate k. The collective Rabi split- 
ting was found to scale as y/~N. All these systems are de- 
scribed by the Dicke model[4( EqnQ]where a single mode 
of photons coupled to an assembly of N atoms with the 
same coupling strength g. 

The importance of various kinds of Dicke models in 
quantum optics ranks the same as the boson Hubbard 
model, Fermionic Hubbard model, Heisenberg model in 
strongly correlated systems and the Ising model in Statis- 
tical mechanics. Since the Dicke model was proposed in 
1954, it was solved in the thermodynamic limit N = 00 
by various methods 0-0] • It was found that when the 
collective atom-photon coupling strength is sufficiently 
large ( Fig.l ), the system gets into a new phase called 
super-radiant phase where there are large number of in- 
verted atoms and also large number photons in the sys- 
tem's ground state [11| . However, so far, there are only 
a few very preliminary exact diagonization (ED) study 
on Dicke models at finite N 0, , its underlying physics 
remains unexplored 10]. It is known that any real sym- 
metry breaking happens only at the thermodynamic limit 
N — > 00, so in principle, there is no real symmetry break- 
ing, so no real super-radiant phase at any finite N. But 
there is a very important new physics for a finite system 
N called quantum phase diffusion in imaginary time at 



finite TV for a continuous symmetry breaking ground state 
at N = 00. The quantum phase diffusion process in a fi- 
nite system is as fundamental and universal as symmetry 
breaking in an infinite system. Here, we will explore the 
quantum phase diffusion process of the Dicke model by 
a 1/7V expansion. We determine the ground state and 
single photon excitation spectrum in both normal and 
superradiant phase. In the normal state, we find a y/~N 
behavior of the collective Rabi splitting in the single pho- 
ton excitation spectrum consistent with the experimen- 
tal data and also determine the corresponding spectral 
weights. In the superradiant phase, we identify a Berry 
phase term which has dramatic effects on both ground 
state and the excitation spectra. The single photon ex- 
citation spectrum has a very low energy quantum phase 
diffusion mode Ed with a high spectral weight and also a 
high energy optical mode E Q with a low spectral weight. 
Their energies and the corresponding spectral weights are 
calculated. The photons are in a number squeezed state. 
The squeezing parameter ( namely, the Mandel Qm fac- 
tor ) is determined. It is the Berry phase which leads 
to the "Sidney Opera " shape in the single photon exci- 
tation spectrum and the consecutive plateaus in photon 
numbers in Fig.l. The Berry phase is also vital to make 
quantitative comparisons between the analytical results 
in this paper and the very preliminary ED results in Q 
and much more extensive ED in Being very strong 
in intensity and has much enhanced signal/noise ratio, 
the number squeezed state from the superradiant phase 
may have wide applications in quantum information pro- 
cessing [l3| and also in the field of high resolution and 
high sensitive measurement [14]. Several experimental 
schemes to realize the superradiant phase of the U(l) 
Dicke model briefly discussed. 

In the U(l) Dicke model a single mode of photons 



2 



Quantum Fluctuation 




FIG. 1: ( Color online ) The single particle excitation 
spectrum E ( in energy unit ) of the U(l) Dicke ( Tavis- 
Cummings ) model Eqn[T] versus the collective atom-photon 
coupling strength g ( in energy unit ). The critical coupling 
is g c — ^/uj c ui a - In the normal state g < g c , the single photon 
spectrum S(lu) has two peaks at the lower branch E- and the 
upper branch E+ with spectral weights c_ ~ 1 and c+ ~ 1 re- 
spectively. The collective Rabi splitting R = E+ — E- ~ g\f~N 
was measured in In the super-radiant phase g > g c , the 
zero mode indicated by the dashed red line at g > g c is lifted 
to the phase diffusion mode Ed in Eqn llOl bv the 1/N correc- 
tion. The single photon spectrum S(ui) has a peak at the low 
frequency phase diffusion mode Ed ~ uj c /N(1 + a) with the 
spectral weight cd ~ 7VA 2 ~ N and a high frequency optical 
mode El = E + Ed with a spectral weight c ~ 1. It is the 
Berry phase ( — 1/2 < a < 1/2 ) effects which lead to the " 
Sidney Opera " shape of the single photon energy spectrum. 



couple to N two level atoms with the same coupling con- 
stant g. The two level atoms can be expressed in terms 
of 3 Pauli matrices a a ,a = 1,2,3. Under the Rotating 
Wave (RW) approximation, the U(l) Dicke model can be 
written as: 
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H u{1) = u, c a!a + -f < + -7= lJ(oVr + h.c.) (1) 

i— 1 v i—l 

where the lu c , uj a are the cavity photon frequency, the en- 
ergy difference of the two atomic levels respectively, the 
g = VNg is the collective photon-atom coupling ( g is 
the individual photon- atom coupling ), the cavity mode 
a could be any one of the two orthogonal polarizations of 
TEMqo cavity modes in [l| . One can also add the atom- 
atom interaction H at - a t to the Eqn[TJ Because H at - a t 
does not change the symmetry of the model, we expect 
all the results achieved in the paper remain qualitatively 
valid. The Hamiltonian Eqn[T] has the U(l) symmetry 
a — » ae l0 ,<T~ — > <j~e 10 . In the normal phase, (a) = 0, 
the U(l) symmetry is respected. In the super-radiant 
phase, (a) ^ 0, the U(l) symmetry is spontaneously bro- 
ken. The model EqnQ]was studied by a large N expan- 
sion in Ref.jH, Q. However, they did not extract any 
important physics at the order of 1/N. In this paper, 



we will show that by carefully analyzing the effects of 
the zero mode in the super-radiant phase, one can ex- 
tract the most important physics of quantum phase dif- 
fusion at the order 1/N. In the large N expansion in the 
magnetic systems [lfj, N is the order of the magnetic 
symmetry group 0(N),SU(N), Sp(2N) with N = 3,2, 1 
respectively. However here N ~ 10 5 is the number of 
atoms, so 1/N ~ 10~ 5 expansion is quite accurate. 

Following standard large N techniques developed in 
fljj . after re-scaling the photon field a — > y/Na, inte- 
grating out the spin degree of freedoms, one can get 
an effective action S e ff[a] in terms of the photon field 
only, then perform a large N expansion. At N — > oo, 
the photon mean field value (a(r)) = A is determined 
0-0| by the saddle point equation: lo c X = g 2 \ tan ^ E 
where E = ^{^f) 2 + g 2 X 2 and f3 = l/k B T is the in- 
verse temperature. At T = 0, in the normal phase 
g < g c — y/uja.ujc, A = 0, the photon number n/N = 0; 
in the super-radiant phase g > g c , A 7^ 0, the photon 
number n/N = X 2 ~ g — g c when g is slightly above g c . 
Its phase diagram and photon number at N — 00 and 
finite N is shown in Fig.l and 2 respectively. 

At finite A, writing a = A + ip where ip describes the 
photon fluctuation around its mean field value A, one can 

expand the S[a] = So + S2 + S3 H to second order [§] 

in 
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where the explicit expressions of the K\ , K2 are given in 
Q, but are not needed in the following. 

In the normal phase g < g c , A = 0. One can see the 
normal Green function at T = in Eqn|2]has two poles 
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T = takes: 



with the spectral weights 
= ( Fig.l ). After the analytic 

it, the one photon Green function at 
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where we also put back the re-scaling factor of the photon 
field a — > \N a. It leads to the two peaks at the two poles 
E± in the single photon energy spectrum shown in the 
Fig.l. At the resonance A c = uj c — uj a = 0, the collective 
Rabi splitting E + — E- = 4g = 4^/Ng ~ y/N shown in 
Fig.l was measured in [l[. Note that due to g a+ > g a _ 
in[l|, the coefficients of the \f~N are different for the two 
different polarizations. The intensity ratio of the two 
peaks — = E+ Z^ a seems has not been measured yet in 

a. c_ " 

However, at the super-radiant phase g > g c , A 7^ 0. the 
anomalous term K2 7^ and |A"i| 2 — |i^2 1 2 contains a zero 
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mode shown as a red dashed line in Fig.l, in addition to 
the pole at a high frequency E Q — \J (oj c + Lo a ) 2 + 4g 2 \ 2 
shown as the blue dashed line in Fig.l. This "zero " mode 
is nothing but the " Goldstone " mode due to the global 
U(l) symmetry breaking in the super-radiant phase. The 
important physics behind this "zero" mode was never 
addressed in the previous literatures [1, [|[ . Here we will 
explore the remarkable properties of this "zero" mode. 
Because of the infra-red divergences from this zero mode, 
the 1 /N expansion in the Cartesian coordinates need to 
be summed to infinite orders to lead to a finite physical 
result. It turns out that the non-perturbative effects of 
the zero mode can be more easily analyzed in the polar 
coordinate ( or phase representation ) by writing a = 
\+4>i+iip2 — \/}^+5pe l6 , then to linear order in Sp and 
9: tp2 — A0. In Eqn|2l by integrating out the massive ipi 
mode, using ip2 = A0, also paying a special attention to 
the Berry phase term coming from the angle variable 
9, one can show that the dynamics of the phase 9 is given 
by: 



S 2 \9} = iNX 2 d T 9 
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2X 2 lu 2 (lu 2 + E 2 ) 
2/3^ w c (w 2 + VA 2 ) 



E 



|0MI 2 (4) 



In the following, we will discuss the the zero mode and 
the optical mode respectively. 

In the low frequency w < £„ limit where the magni- 
tude fluctuations can be dropped, Eqn|4] reduces to: 

C PD [9]=iN\ 2 d T 9+^{d T 9) 2 = ^L(d T 9 + taD) 2 (5) 

with the quantum phase diffusion constant D = E i-L . 

In the EqnfSl we have denoted TV A 2 = No+a where No — 
\N\ 2 \ is the closest integer to iVA 2 , so -1/2 < a < 1/2. 

The corresponding quantum phase diffusion Hamilto- 
nian is: 



Hpd[9] = j(6N ph -a) 2 



(6) 



where SN p h — N p h — Nq is the photon number fluctuation 
around its ground state value Nq and is conjugate to the 
phase 9: [9,SN p h] — ih. In fact, Eqnj6]can be considered 
as the Hamiltonian of a particle moving along a ring with 
a very large inertial of moment I = 1/D subject to a 
fractional flux / = (f>/(j)o = ct. 

In Eqn[5l after defining 9(t) = 9(t) + iaDr, one can 
easily show that 
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which is a phase diffusion in imaginary time t 12|, [18 1 
with the phase diffusion constant D. Only in the thermo- 
dynamic limit N — > oo, a state with a given initial phase 
will stick to this phase as the time evolves, so we have a 
spontaneously broken U(l) symmetry. However, for any 
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FIG. 2: ( Color online ) The average photon number in the 
ground state at N — oo and at a finite N < oo versus the 
collective atom-photon coupling strength g ( in energy unit ). 
The photon is very close to be the Fock state as shown by the 
steps in the superradiant phase at N < oo. It is the Berry 
phase (— l/2<a<l/2) effects which lead to the steps in 
the photon number. 



finite N, the initial phase has to diffuse with the phase 
diffusion constant D ~ 1/N. The diffusion time scale in 
the imaginary time beyond which there is no more phase 
coherence is td = 1 / D ~ N/uj c which is finite for any fi- 
nite N. This can also be called phase "de-coherence" 
time in the imaginary time [l8j ]. 

From EqnfSl it is easy to see that the gapless nature of 
the phase diffusion mode in EqnjS] leads to the vanishing 
of the order parameter inside the super-radiant phase: 



(a) =0 



(8) 



So the U(l) symmetry is restored by the phase diffusion. 
From Eqnf5l after doing the analytic continuation t — > it, 
we can get: 



(a t (t)a(0)) s = N\ 2 e-^ +a ^ Dt 



(9) 



where we also put back the re-scaling factor of the photon 
field a — > \/Na. It is also easy to see that (a(t)a(0))s = 0, 
so there is no quadrature squeezing anymore at any finite 
N [l9j . In fact, all these results can also be achieved by 
using the Hamiltonian EqnjBJ 

Eqn[9]leads to the result that the energy of the "zero 
energy mode " ( Goldstone mode ) at N — oo was " lifted" 
to a quantum phase "diffusion " mode at any finite N 
with a finite small positive frequency fl8"i ]: 



e d = (r+aO-D = 77 — — — r- 1 

2 [[U3 C + w 0/ 



Ag 2 X 2 ]N 



ujJN (10) 



It is the Berry phase effect which leads to the periodic 
jumps in the Fig.l. The Fourier transform of Eqn|H]leads 
to the Fluorescence spectrum S{ui) = N\ 2 8(lo—Ed) with 
the spectral weight c° D = NX 2 ~ N. 

Now we study the photon statistics. If neglecting 
the magnitude fluctuation, the quantum phase diffusion 
Hamiltonian Eqn[6] shows that the ground state is a pho- 
ton Fock state with eigenvalue Nq which jumps by 1 in 
all the plateaus ending at a = 1/2 in the Fig. 2. Now 
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we incorporate the magnitude fluctuation. In the Eqn|2] 
by integrating out the imaginary part ip2(u) and using 
ipi = Sp/2X, one can get the effective action for the mag- 
nitude fluctuations 



C 2 (5p) = N- 



E 



8X 2 u, 



(ii) 



where we find the Mandel factor Qm 



{(8Nf)-(N) _ 
(N) 

-g- so the deviation from the Fock state at any 
given plateau in Fig. 2 is given by Because =^f- <C I 
in the g ^> g c limit, it is very close to be a Fock state. 
This is a highly non-classical state with Sub-Poissonian 
photon statistics. It has very strong signal (N p h) = NX 2 , 
but nearly no photon number noise, so it has a very large 
signal to noise ratio which could be crucial for quantum 
information processing [13J and also in the field of high 
resolution and high sensitive measurement 14j. 



Using the polar representation a — \/X 2 + Spe 10 ~ 
Xe l9 + ^-e l9 + 0{l/N), one can evaluate the photon cor- 
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relation function: 



(Ta t (r)a(0)) =N\ 2 (e 
N 



-i(9(r)-e( 



0))\ 



+ —(Sp(r)Sp(0))(e-^- s ^)+O(l/N)(l2) 

where the T means imaginary time ordered. By eval- 
uating the first and second term from Eqns|3] and II H 
we can identify not only the quantum phase diffusion 
mode Ed = D{\ + a) with the corresponding spec- 
tral weight c D — NX 2 — ^^rj^r^ , but also the opti- 
cal mode £"* = E a + Ed with the corresponding spec- 
tral weight c = ^(^+ 3 "°) 2 + 2j£ ( Fig.l ). Note that 

E Q — Ed = E a in Fig.l is independent of the Berry phase 
a. So the total energy in the optical frequency peak 
~ E Q x c ~ u c is comparable to that in the phase diffu- 
sion mode ~ Ed x cd ~ u c /N x iV~ ui c . 

If one introduce the total "spin" of the N two level 

atoms J z = J2i °f> J+ = Hi a t> J ~ = Hi 17 i and con- 
fine the Hilbcrt space only to J = N/2, then EqnfT]can 
be simplified to the J—U(l) Dicke model which was stud- 
ied by an exact diagonization (ED) in Q. The authors 
in Q found that there are a series of ground state en- 
ergy level crossings as the g gets into the super-radiant 
regime and interpreted them as consecutive "quantum 
phase transitions " . But they did not study any excited 
states. In (l6| . we performed a much more extensive ED 
study directly on the U(l) Dicke model EqnfT] and not 
only calculated the ground states, but also all the ex- 
cited energy levels. We also identified a series of ground 
state energy level crossings in the super-radiant regime. 
By comparing with the analytic results achieved in this 
paper, we found all these ground states crossings are pre- 
cisely due to the periodic changes of the Berry phase a in 
Eqns l4"I5Io1 They are not consecutive "quantum phase- 
like transitions " as claimed in [7fl . We also found one to 



one quantitative matches between the low energy phase 
diffusion mode Ed, also the high energy optical mode 
E Q = E + E D and the excited levels found by the ED 
in [16| at iV as small N = 5. The complete comparisons 
will be presented in [lj|. 

It remains experimentally challenging to move into the 
superradiant regime which requires the collective photon- 
atom coupling g — \Z~Ng > g c = \fu c w a ~ 2n x 10 5 GHz 
for the optical cavity used in [1(. The collective Rabi 
splitting y/Ng ~ 20GHz in [l[ is still much smaller than 
= 2tt x 10 5 GHz, so not even close to the superradi 



ant regime in Fig.l. It was proposed in 2(3- 22 1 that the 
super-radiant regime can be realized by using a cavity- 
plus-laser-mediated Raman transitions between a pair of 
stable atomic ground states, therefore also suppress the 
spontaneous emission 7. All the parameters in EqnJT]can 
be controlled by the external laser frequencies and inten- 
sities, so the characteristics energy scales in the effective 
two level atoms are no longer those of optical photons and 
dipole coupling, but those associated with Raman tran- 
sition rates and light shifts. Indeed, using this scheme, 
the super-radiant phase in the Z 2 Dicke model [j| |20| 
was reached by using both thermal atoms 2l| and the 
cold atoms in the BEC (22|. We expect this scheme 
may also be used to realize the super-radiant phase of 
the U(l) Dicke model shown in Fig.l. Because the mi- 
crowave circuit cavity has much lower cavity frequency 
and the individual photon-qubit g can also be made very 
large, so the superradiant phase could also be realized in 
superconducting qubits or quantum dots inside a circuit 
cavity in the future. In the experiments, there is also a 
weak dissipation n <§; g. In a future publication, follow- 
ing the procedures in [19(, we will study the effects of k 
on the number squeezed state. 
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